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Recent measurements of the force versus extension curves in 
10 ' stretched single stranded DNA, under conditions where the hy- 
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drogen bonding between complementary bases is inhibited, pro- 
vide a new handle for the study of self-avoiding effects in a flex- 
es \ ible polymer. We report in this paper upon analytic compu- 
tations of the force versus extension curves within a continu- 
ous version of the freely joining chain model, with monomer- 
monomer repulsive interactions. The problem is formulated as a 
S ■ Statistical Field Theory model, endowed with a fixed cutoff as- 
. sociated with the curvature of the extension versus force curve, 
P3 \ in the free polymer limit. Using a Field Theory version of the 
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Hartree-Fock approximation, the self-avoiding single polymer 
problem reduces to the solving of a one-dimension self-consistent 
integral equation. Taking the short range potential limit, we 
^ ' obtain an explicit analytical solution, involving the roots of an 

■ algebraic equation. In the low force regime, we find that the 

slope of the relative extension at zero force increases steadily 
with the total monomer number N , while it stays constant for 
a free polymer. We compare, for = 100, our analytic results 
with those obtained by Monte Carlo simulations of a discrete 
freely joining chain, in presence of excluded volume effects. De- 
spite a difference between the repulsive potential ranges d, an 
appropriate choice of our single free parameter leads to a good 
agreement between the two computations. This may be consid- 
ered as an indication that the shape of the extension curve is 
not very sensitive to d, if it is equal to or less than the monomer 
length. 
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Introduction and Summary 

During tlie last few years, ttie elastic properties of single biopolymer molecules, un- 
der physical conditions close to those encountered in living organisms, have been the 
subject of extensive investigations |l|, 0, H- In particular, the force versus extension 
curves of double stranded DNA (dsDNA) have been measured accurately within a 
wide range of pulling forces. In absence of supercoiling constraints, the entropic elas- 
ticity of the dsDNA molecule is very well described 0] by an elastic string Q model, 
the so-called Worm Like Chain (WLC) model, involving a single rigidity associated 
with the resistance against the polymer bending. The dsDNA persistence length A is 
given, within the WLC model, by the bending rigidity written with appropriate units. 
Micromanipulations experiments lead to values clustering around A ~ 50 nm, which 
corresponds to 150 base pairs, hence the qualification "stiff' polymer sometimes given 
to the dsDNA. The monomer-monomer interactions in dsDNA are dominated by the 
screened Coulomb repulsion and they do not seem to lead to significant corrections 
to the WLC model predictions]^], with the exception of very long chains, longer than 
50 /im§. 

There has been recently a growing interest for micromanipulation experiments 



1^, |I0| involving the stretching of single stranded (ss)DNA. The elasticities of ssDNA 
and dsDNA are expected to be different, for at least two reasons. One, despite 
the lack of any accurate determination, it is commonly believed that the ssDNA 
persistence lies in the range 0.7 to 3 nmP], corresponding roughly to a base pair 
number between 2 and 7. In presence of thermal fluctuations, the molecular local 
axis loses the memory of its direction after few base pairs. We express this fact by 
saying that the ssDNA is a "locally flexible" polymer". In such a case, the WLC 
model does not lead to an adequate description of the "free" ssDNA molecular chain. 
A more realistic picture is provided by the freely joining chain (FJC) model, where 
the relative direction of two adjacent elementary links is allowed to fluctuate freely. 
As a consequence of the ssDNA high flexibility, the self avoiding-effects are expected 
to be more important than in the case of a "stiff" polymer like dsDNA, if the physical 
environment is such that the range of the repulsive monomer-monomer interaction is 
about few nm. The second reason for a difference between the ssDNA and dsDNA 
entropic elasticities is obviously the attractive interaction between complementary 
bases, leading to the so-called "hairpin" structures which have to be opened up in 
order to stretch the molecular chain. A theoretical analysis of this mechanism has 



been performed recently |TT|, with the help of the "rainbow" approximation. The 



results are in good agreement with the observations of ref. [|T^ when the buffer solution 
contains cations which stabilize the hairpin structures. 



2 



Conversely, ssDNA stretching experiments |jTO[ have been performed under phys- 
ical conditions where hydrogen bonding interactions between complementary bases 
are suppressed, leaving the screened Coulomb repulsion as the dominant monomer- 
monomer interaction. The necessity of accounting properly for self-avoiding effects 
appears clearly from the lack of agreement between the experimental data and the 
ideal polymer model predictions, like the FJC model. This is confirmed [|10| by a 
Monte-Carlo (MC) simulation of the force vs extension curves within a FJC model 
with excluded volume effects. 

In this paper, we compute, by Statistical Field Theory methods, the entropic 
elasticity of locally flexible polymers in presence of repulsive monomer-monomer in- 
teractions. 

As a first step (section 1), we have constructed a continuous version of the FJC 
model, taking the monomer number n as the running variable along the chain. It 
is well known that one can map the free polymer problem upon a single particle 
quantum problem in an Euclidian space-time. In the present case, the imaginary 
time is just —i x n. The partition function matrix , relative to arbitrary states of 
the molecule free ends, is given hy Z = exp — Hq , where Hq is the Hamiltonian 
of the attached quantum problem and A^ the total monomer number. In the case of 
the continuous version of FJC model, the Hamiltonian is purely kinetic in the zero 

force limit: Hq = ho(^p^, where p = — z V is the momentum operator conjugated 

to the vector f = — tq, joining the two ends of the chain. The force stretching 
energy —f- F is then implemented through the replacement: p — > p — iF. One proves 
easily the following simple relation between the relative extension vs force function 
Co{F) and the fictitious free particle Hamiltonian: Co(-^) = 2-F/io(— -F^). In principle, 
the Hamiltonian Hq is then obtained by integrating the extension vs force function 
along a finite interval along the imaginary F axis. Unfortunately, the things are 
not so simple since Cfjc{F) has poles along integration path. We have been able 
to bypass this difficulty by building a function Cr{F) regular along the imaginary F 
axis and differing from Cfjc{F) by less than a few % along the whole real F axis. 
The coincidence is even better than 1% in the low to medium force regime where the 
self- avoiding effects are the most important. The Hamiltonian Hq is written finally 
as Hq = Hq (^{p — iPy^ where and fiQip^) is an analytic function of p increasing faster 

than p^. In our construction, the quartic expansion h^\p'^) = + is in fact 
obtained by integrating the cubic expansion of (pjciF): where b is the length of the 
FJC elementary link. The quadratic term corresponds to the Gaussian model; when 
it is used in association with a short range (Dirac 5 function) momomer-monomer 
repulsive interaction, it is called the Edwards model. The quartic terms provides a 
natural cutoff A = 306~^ for the divergences appearing in the perturbation expansion 
of the Edwards model with respect to the monomer-monomer interactions. It will 
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turn out that in the exploration of the low to medium force range < {3 F b < 1.5 
the quartic Hamiltonian hl^\p'^) can be substituted to h^lp'^), allowing for analytic 
computations. 

In section 2 we develop an Hartree-Fock treatment of self-avoiding effects within 
our continuous version of the FJC model. With the help of a well known functional 
integral identity ||13||, the self-avoiding polymer problem is exactly mapped into that of 
a quantum single particle moving in a stochastic imaginary potential, described by the 
Hamiltonian: = Hq+z g 0(r). The functional probability measure of the field 0(r) 
is of the Gaussian type, specified by the two fields product average: ( (f>{fl) 0(^1) )^ = 
V{\fi — r^l) where V{\fl — "t^D is the repulsive monomer-monomer potential. The 
partition function of the self-avoiding polymer with fixed free ends is given by the 
stochastic field average: Z{fN — tq, F, N) = (r^r | (exp — (^Hq + i g )^\r\)). To 

proceed, it is convenient to perform the Laplace transform of the partition function at 
zero force Z{f,0,N) z{f,T), where r is the variable conjugated to A^. The result 

takes the following simple form, z{f, r) = (r | ( ^/iq ^ + ig + ''"^ ) I 0)) which 

can be readily written as a Feynman diagram expansion in power of g"^. It allows us 
to use some of the Quantum Field Theory machinery, in particular the Dyson integral 
equations. Following a standard procedure |jl2|, we introduce an approximate kernel 
in the Dyson integral equation for the propagator. We arrive in a rather direct way 
to the Hartree-Fock equation which performs the exact summation of the "rainbow" 
diagrams. The self-consistent integral equation takes a remarkable simple form in the 
Fourier space. It involves then the Fourier transform z{p, r) of 2;(r, r), where p is the 
momentum conjugated to the relative free end coordinate f. 

^ ^ [d'qz(q,r)V{\p-q\) (1) 



z{p,t) Zo{p,t) 2tt 

where V"(|g|) is the Fourier transform of the potential V^(|r|) and zq{p, t) corresponds 
to the free polymer partition function. 

In this exploratory paper we have used a zero range potential, proportional to a 
Dirac 5 function. The Hartree-Fock integral equation reduces in this limit to a nu- 
merical equation. We will show that z{p, r) can be written as (/iQ(p^)+/i(r))"-'^, where 
the function /i(r) obeys the self-consistent equation: /i(r) = r + ^F{k, 6yu(r)), where 
K = P and A a positive dimensionless constant measuring the strength of the re- 
pulsive monomer-monomer potential. In the low to medium force regime, F{k, t) is a 
rather simple algebraic function of r. The extension versus vs force function ({F, N) 
for the self avoiding polymer is then obtained by taking the logarithmic derivative 
of the inverse Laplace transform of (/io(— F^) + /i(r))~^. The involved integral is 
dominated by the residues of the poles at the roots of /iq(— F^) + yu(r) = 0. 
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In the low force regime < (3 F b < 0.5, the computation of the extension vs 
force function ({F,N) can be performed analytically up to the very end. The final 
expression is not by itself very illuminating, so we have plotted a set of extension 
curves for a wide range of monomer number: 50 < iV < 5000. The most remarkable 
feature is the steady increase with of the relative extension vs force slope at the 
origin {F = 0). This contrasts with the case of a free polymer where this quantity 
is independent of N. This behaviour is also predicted in the Edwards model using 
scaling arguments ||16||, but the effect is notably smaller than in our approach. Such a 
difference was to be expected since Renormalization Group methods cannot be applied 
to our model which is endowed with a fixed cut-off, associated with the extension vs 
force curvature of the free polymer. 

The section 4 is devoted to a comparison of the analytical results of our Hartree- 
Fock method with those obtained by a MC simulation |jlO| of a discrete freely joining 



chain, with excluded volume effects. In both cases, the total number N of monomer is 
one hundred. Our analytic treatment is strictly valid in the low to medium force range 

< (3 F b < 1.5, but we have extended our results to higher forces with the help of a 
reasonably safe extrapolation procedure. We are able in this way to reach the force 
domain, (3Fb > 3.5, where self-avoiding effects are disappearing. An appropriate 
choice of our single adjustable parameter, A = 0.7, leads to a predicted extension vs 
force curve which is in remarkable agreement with the MC simulations. It should 
be stressed that there is a notable difference between the physical inputs in the two 
approaches. In our exploratory computation, we use a zero range potential while 
the MC simulations were performed with an excluded volume radius d equal to the 
monomer length b. This may suggest that the shape of the extension vs force curve 
is not very sensitive the value d of the repulsive potential range ii d <b. 

1 A continuous formalism for a free flexible polymer. 

In this paper, the partition function for a stretched " locally flexible polymer" will 
be given, in absence of self-avoiding effects, by the following path integral: 

N 



Zo{fN, ro, F) = J V[r\ exp - ^ dn So{n 



where n stands for the monomer number which is treated as a continuous variable 
as a result of a coarse-graining of the polymer chain ; Soin) is the elastic energy per 
polymer written in thermal energy unit fc^T.The function eo(a;) is assumed here to 
be analytic. The equation could be considered as the definition of what we mean 
by a " locally flexible polymer" . In the literature about polymer self-avoiding effects. 
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the function eo(a;) is often taken to be linear: £o{n) = ^{^Y ~ ^ ' ^ where b is the 
root mean square distance of two adjacent polymers. The difficulty with this model- 
known as the "Gaussian model"- lies in the fact that the extension vs force curve 
turns out to be a straight line, which is quite unrealistic. The philosophy taken up 
in this paper is to consider as the true starting point the extension vs force function 
in absence of self-avoidance : < zn > /N = Co{F)i rather than the elastic energy 
eo ((^)^) appearing in equation d). 

1.1 An Hamiltonian formalism to relate the extension vs force function to the linear 
elasticity density. 

We are now going to write the discretized form of the partition function given in 
equation (|^) in terms of the transfer matrix operator Tq : 

- N 

Zo{rN,ro) = (rVlTo \ro), 

{rn+l\fo\fn) = exp - (eo {^{fn+l - fnf) - F ■ {fn+1 - fn)). (3) 

The operator Tq is invariant upon space translations and is then expected to be 

3 

diagonal within the plane wave basis: {rip) = {2tt)~^ exp{ip ■ r). The transfer 
matrix in the momentum space is then readily obtained : 

(Pn+l| To \Pn) = S^{Pn+l " Pn) Tq - iF)^) , (4) 

where Tq {p^) is given in terms of eo (r^) by the following Fourier integral: 

To (p^^ = J (i^rexpzp ■ r exp eo(r^)) . (5) 

Introducing the momentum operator p = —iV the transfer matrix operator can be 
written under the simple form: Tq = Tq (j^p — iPy^. We arrive in this way to a 

compact formula for the partition function, which will lead us to the Hamiltonian Hq 
of an associated quantum mechanical problem: 

ZoifN - ro, F, N) = (fjvl {ip- iFf) |Fo) = (riv| exp(-OTo) |ro). 

Ho = ho [(p - iFf)= -log To [(p-iFf), (6) 

where in our writing of the partition function we have made explicit the fact that 
is invariant under space translation. In the simple case of the Gaussian model, one 
finds: ^0 = -Y (V + Ff. 
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Let us proceed to the evaluation of the relative extension vs force function Co(-^)- 
The extension average < r > is given by : 

<f>= Aiog(/ dr^Z,{r,F,N)] = -N^h,{-F^). (7) 
oF \J / oF 

Taking the stretching force along the z axis, one gets the relation which will allow us 
to construct the effective Hamiltonian Hq from the extension vs force curve Co(/)- 

UF) = — = 2Fh',{-F') = -2f- 



N ' 'roi-F^y 

Let us now consider the case of the freely joining chain (FJC) which could be con- 
sidered as a reasonable starting point for the study of self-avoiding effects for a " 
locally flexible" polymer. The transfer matrix operator for a chain of freely join- 
ing monomers, having a fixed length b, can be written as: {rn+i\To{FJC) |r„) = 
6 — fn\ — b) exp F ■ f. It can be easily written in momentum space form as in 

eq. d^) by taking: Tq (p^) = sm{pb)/{pb). Inserting this result in eq. (H), one gets 
immediately the FJC model extension vs force curve in terms of the reduced force 
f = Fb: 

CfjcU) = b (coth(/) -j]={-^ + 0{f). (9) 



1.2 Construction of an effective regular Hamiltonian from the extension vs force 
function. 

If one applies blindly the formula (^) in order to get the effective Hamiltonian Ho, 
one gets /io(p^) = — log C^™^^ ), which is clearly singular when pb = nn, n being 
an arbitrary non zero integer. One may try to regularize the model by replacing 
in (fVi+i| Tq{FJC) \ fn) the 6 function by a Gaussian, allowing for fluctuations of the 
monomer length b. This does not solve really the problem since To(p^) still oscillates 
around zero and holp"^) does not increase fast enough with p^. The fact that the above 
formal derivation of hoijp) leads to a singular result should not be considered as a 
surprise. A typical raw FJC contour line, in absence of any coarse-graining, is very 
far from what could be obtained by the discretization of a 3D-curve with smoothly 
varying derivatives of the running coordinate ^. 

We have adopted here the view that there is nothing sacred about the FJC model: 
it is just a guide to get a reasonable extension vs force curve for a "locally flexible" 
polymer. 

Let us sketch the construction of a regular Hamiltonian h^ip^) associated with 
an extension vs force function Cr(/) differing from (pjcif) by less than few % for 
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Extens ion 




Figure 1: This figure gives the results of a fit of the extension vs force function (pjcif) by a 
superposition of two error functions. The idea behind this procedure is to build a continuous variant 
of the FJC mode, where the associated free Hamiltonian is regularized, while keeping the essential 
features of the FJC model elasticity. We see clearly on this figure that this goal has been successfully 
achieved. We have required the fitted function Cr(/), with the label " Fit" to coincide with CfjcH) 
to third order in the reduced force f = F/b. The dotted curve, with the label "Cubic", represents 
the third order polynomial expansion of CfjcH)- 

the whole range of values of the reduced force /. The result is exhibited on Fig.|l[ 
A convenient starting point is the derivative Cpjcif) ~ ~sinh~^(/) + /~^. It is a 
bell-shaped even function of / we write as x{P)- the case of the FJC model 
x{z) is an analytic function of z with an infinite set of poles at Zn = — (nvr)^, which 
are responsible for the singularities of holp"^). We shall suppress them by replacing 
Xfjc{z) by a function regular along the whole z, having the general form: Xr{z) = 

dm g{m) exp{—mz). Here g{m) is a distribution which is non zero if and only 
if the real number m belongs to the finite domain < m < M where M is a finite 
positive real. Note that this representation is not valid for Xfjc{z) unless M is infinite. 

As in all regularization methods, there is always some arbitrariness. Two cri- 
terion's have to be used as guiding principles: one is mathematical simplicity, the 
second is keeping to a minimum the number of ad hoc parameters. We have found 
that a very simple form of the distribution g{m) is adequate for our purpose: Xr{z) 
is taken as the sum of two exponentials. Returning to the / variable, it implies that 
the regularized bell-shaped curve C'(/) is the sum of two Gaussian. The main effect 
of the regularization is to change the shape of the high force tail of C'(/) ^ domain 
where self-avoiding effects are becoming very weak. One of the peculiar feature of 
the analytic continuation f ^ ip is that a modification performed for large positive 
/ values can induce or suppress singularities for finite values of p : vr, 27r 
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By performing a simple quadrature, we arrive finally at the following expression 
for the regularized extension function Cr(/) : 

Cr(/) = ^ ((1 - c) erf (ai /) + (! + c) erf (a^ /)) , (10) 

where erf (2;) = /q°° dt exp —t"^. Note that if ai and 02 are positive, then Cr(oo) = 1. 
The determination of the parameters ai, a2 and c is performed in two steps: 

Step 1. ai and 02 are obtained as functions of c by requiring the equality of 
the third order expansions of Cr(/) and (pjcif) with respect to /. As it is clear from 
Fig.|l], this constraint insures that Cr(/) and (pjcif) coincide in the low to medium 
force range < / < 1.5, where self- avoiding effects are expected to be important. 

Step 2. The parameter c results from a minimization of the mean square aver- 
age difference : A2(/™) = j;: St df {Uf) ' CfjcU)?- For |/| > 10 the molecule 
extension is above 90% of its maximum. Self-avoiding effects are thus expected to be 
small, and this is confirmed by a look at the curves displayed on FIG. |^. This justifies 
our choice: fm = 10. The resulting numerical values for the constrained parameters 
are then given by : 

c = -0.035 , ai = 0.093 , aa = 0.484 (11) 

With the above numbers, one gets for the root mean square difference : ^A2{fm) = 
0.0198. Such a small difference gives a measure of the success of our enterprise to 
suppress the singularities of /io(p^) while keeping to a few percent level the modifica- 
tions of the FJC model extension vs force function. It is also important to note that 
there is no free parameter left. 

We are now ready to write down the effective Hamiltonian Hq which will be used 
in this work : 

Ho 



where erfi(z) = eTf{iz)/i. The function h^ip^) is displayed on Fig.^ 

We have verified that the use of h^ip^) does not introduce any significant non- 
physical artefact in the model. To do that, we have worked backward: taking h^ip^) 
as starting point, we have computed the probability distribution of /„ = l^n+i — rn\, 
given by P(/„) = ll f(/„) = 1^ exp -eo(/^). Using eq.(D, we get f(/„) by performing 
the 3Z) inverse Fourier transform of exp— /iq(p^). Discarding a very weak damped 
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hl[(p-tF) 
CP ^ 

du (r {i u 



-(1 — c) ( — L(l — exp(aip^6^) +p6erfi(aip6) 



+ (c ^ -c, ai — > 02) 



(12) 



Hamiltonian 




Figure 2: On this figure, the solid line represents the effective Hamiltonian ft,Q(p^), obtained by a 
simple quadrature from the analytic continuation / — > of the force -extension curve Crif)- The 
dashed curve gives the fourth order expansion: h'^\p'^) = ^-^ — h ^j^- It nearly coincides with 
hgip^) in the interval < p < 3 and will be used to get the analytic results given in this paper. Note 
the fast increase of two the curves for large values of . It guarantees the absence of divergences in 
the field theory description of self-avoiding effects which will be adopted in the present work. 



oscillating tail at large /„, having a maximum amplitude which stays below the 1.5% 
level, we get a well behaved Gaussian-like bell- shaped curve centered at /„ ~ 6 with 
a half-width A/„ ^ 6/2. In our procedure, regularity conditions are enforced at each 
step. Our derived probability distribution P(/„) is not expected to be a Dirac 6 
function but should rather be a peaked curve with a finite width, which should be 
of the order of the coarse-graining length resolution needed to wash out the strong 
fluctuations of the raw FJC contour line. 

The quartic expansion of /lo(p^) is, by construction, independent of the parameters 
Oi, a2 and c : 

.?V)^!^.!^ (13) 

As it is apparent on Fig ^ and 0, it could constitute an adequate starting point for 
the study of the self-avoiding effects of the FJC model in the force range < / < 1.5, 
which will be studied in the next sections. The effective Hamiltonian Hq^^ associated 
with hf^\p'^) can be viewed as a cutoff version of the Gaussian model. Indeed, the 
corresponding Green function in the momentum space: Go{p'^,z) = {z — hQ'\p'^))~^ 
behaves as when S> = 30 h~^. The cut-off A is not here an ad hoc ingredient, 
since it is related to a simple physical quantity: the curvature of the extension vs force 
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curve in the low force regime. 

We would like to stress that h^o\—p) is not valid beyond / > 2. Higher poly- 
nomial expansion of /io(p^) = — log(sin(p)/p), as a substitute for /io(p^), will lead to 
difficulties due to the divergence of the entire series expansion of /io(p^) a.t p = n, 2n... 



2 The Hartree-Fock approximation for Self-avoiding Polymers. 

Let us incorporate in the partition function the self-avoiding effects induced by the 
monomer-monomer repulsive central potential g"^ V^d^i — ^'2!) '■ 

Z{fN-fo,F,N) = Jv[r\exp-J^'' dn (^So{n) + £ dn' V{\f,, - r,,,\)y (14) 

Before entering into the details of our approximation scheme, it useful to derive a 
convenient formula for the evaluation of the chain extension vs force. We introduce the 
space average partition function Z{F^, N) obtained by integrating upon the distance 
r = tn — "^0 between the two free ends of the molecular chain.: 

Z{F^,N) = J dr^ Z{f,F,N) = J dr^ exp{f ■ F)Z{f,0, N), (15) 

where have factorized the stretching potential energy and exploited the rotation in- 
variance of Z{f, 0, A^). Taking the stretching force along the z axis, F = z F, we can 
readily obtained the relative extension vs force expression : 

aF,N) = = -^10^ {Z(F\N)). (16) 

The computation method developed in this section will give in a rather direct way 
the Fourier transform of the partition function of the unstretched polymer: 

Z{k'^,N) = J dPexp{if-k) Z{f,0,N). (17) 

It is then of interest to write the extention vs force formula in terms of this quantity. 
To do that, let us take again F along the z axis and perform the analytic continuation 



F ^ ip upon the r.h.s. of eq. (|15|). Inserting the Fourier expansion of Z(r, 0, A^), 
one can readily integrate over r and arrive at the relation Z{p'^,N) = Z{—p'^,N). 
Performing the analjd^ic continuation backwards, one gets the searched-for formula: 

C(^',iv) = ^ = ^log(^(-F^iv 
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2.1 A Statistical Field Theory model for self-avoiding effects in flexible polymers 



Following a standard procedure [13, 14], we are going to transform the evaluation of 
the r.h.s. of (0 ) into a field theory problem involving an auxiliary scalar field 0(r). 
The basic tool is the well known Gaussian functional integral identity : 

I rN rN 



£[V] = exp--/ dn dn' g^V{\fn-rn>\) (19) 
I Jo Jo 

= (det Vy/^J V[<j)]exp- (^J^'' dmg<f){rn)^ 

X exp - (^J d^ri d^ r2 ^V'^in - ^2)0(n) 0(^*2)^ . (20) 

In the case of the Yukawa potential V{r) = ^^exp— r/a, the inverse operator is 
given by : 

V^\n - ra) = 5%n - r2)i-Vi + a~^). 

We arrive in this way to the Field Theory formula for the partition function upon 
which we are going to build our approximation scheme: 

Z{rN-fo,F,N) = lv[(l)]exp-(^l d^r,d^r2^V~\n-f2)(l)in)(P{f2)^ 

x(rW| exp-OT(F,0) |fo). (21) 
The Hamiltonian H{F,(f)) involves the stochastic imaginary potential ig(f){f):, 

H{F,(P) = ho {$-^Ff)+^g4>{f), (22) 

where the "kinetic" term /lo(p^) is to be identified with h'^\p'^) ( eq. (P^), or even- 
tually with /i[;(p2) ( eq. (|12D). 

2.2 A simple Feynman graph expansion for the Laplace transform of the partition 
function. 

We proceed by performing the Laplace transform of (|2T|) with respect to A^: 

z{fN-ro,F,T) = / dN exp{-TN)Z{fN-ro,F,N) 
Jo 

= {fN\{{ho{iP~^FY)+^gm + Ty')^\fo), (23) 

where we have introduced a compact notation to describe the average of a given 
functional J^[4>] over the stochastic field 0(r): 

( -^[0] )4> = 1 1^m<P] exp - (/ d'r, d' lv~\n - r2)0(r-'i) ^(rl)) . (24) 
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One sees immediately that the average product of an odd number of fields 0(ri) is 
zero. We give some simple examples of even products: 



( 0(rl) 0(rl) )o = V{\r-{-r2\) 
{ 0(r"I) 0(r^2) 0(r"5) (pin) )q = Vdn - r^sl) V{\f2 - n\) + - fl\) V{\r2 - rl|) 

+Vi\r-l-r-l\)Vi\r-2-r-l\) (25) 

The rule is easily generalized to a product of an arbitrary even number of fields 
and leads readily to a Feynman graph expansion of z{f, F, r) in power of g^. The 
perturbation expansion of the resolvant R{t) = {Hq + U + r)^^ ( the operator If is 
defined by {r2\U\fl) = ig5^{f{ — r^) 0(^1) ) is given by the geometrical series, keeping 
only the relevant even terms : 

R{t) = Ro{r) + Mt) tj Ro{t) If Ro{t) 

+Roir)U Roir) U i?o(r) U i?o(r) U i?o(r) + (26) 

where Ro{t) = (Hq + t)^^ . To make more apparent the analogy with the Field Theory 
formalism we define the Green functions: 

Go{r2 -ri,T) = Zo{r2 - n, F, r) = (r^ |i?o(^) In) 

G{f2-n,T) = zif2-n,F,T) = {r-2\{RiT))^\n) (27) 

2.3 An Hartree-Fock approximation from the Dyson equation. 

Following Dyson [|l^, we introduce the proper self-energy S(r2 — fi) obtained by 



summing all self-energy diagrams which cannot be divided into two disconnected 
pieces by cutting the "particle line" associated with the Green functions Go(^2— "H,, t). 
Pure topological arguments lead to the well known Dyson equation for the propagator: 

G{f^2-ri,r) = Go(^2-n,i") + 

J £rsd\,Goif2 - rs, r)S(f3 - f^)Gif^ - n, r). (28) 

Taken in their exact form, the Dyson equations have played a very important role 
in the Renormalization Theory but they are of little use for practical computations. 
They are, however, very convenient to devise non trivial infinite summation scheme 
by taking an approximate form of the proper self-energy. The best known example of 
this procedure is the Hartree-Fock (HF) method [| which relies upon the approximate 
formula: 

S(f2 - fi) ^ S«(f2 - fi) = -g^ Gif2 - fi, r) Vif2 - n) (29) 

^ The formalism used here takes its inspiration from the Fetter and Walecka remarkable book 
dealing, with Quantum Theory of Many Particle Systems [^2) . 
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Figure 3: Graphical representation of the Hartree-Fock approximation. The thin hnc and the thick 
hnes stand respectively for the free and the self-avoiding polymer Green functions. The dashed 
half-ellipses represent the self- avoiding interactions. The two upper diagrams schematize the HF 
approximation used in this paper which performs the summation of the planar graphs. An infinite 
set of non-planar diagrams is generated by including the contribution of the bottom diagram. 



A physical interpretation of the HF approximation can be given by noting that 
5](i)(f*2 — Ti) is obtained by replacing, in the lowest order proper self-energy, the free 
polymer propagator GQ{f,T) by G{f,T) which incorporates in a self consistent way 
higher order monomer-monomer interactions . Although technically rather different, 
the present HF method bears a close analogy to the "rainbow" diagrams summa- 
tion method used to compute "hairpin" effects in a single stranded discretized DNA 
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One of the virtues of the Field Theory formulation of the Hartree-Fock method is 
to allow a rather systematic treatment of the correcting terms. The first step consists 
in adding to the HF proper self-energy given by eq.(^) a correction term S*^^^(r2 — ri) 
which will generate an infinite set of non-planar diagrams. It is constructed from the 
lowest order non-planar diagram contribution to the proper self-energy in the same 
way as T,HF{r2 — ri) was obtained from the lowest order proper self energy: 

S^2)(f2 - fi) = J d^r^d^riG{f2 - n)y{r2 - x 

G{n - rs)Vin - ri)G(f3 - n). (30) 

We give in FIG.|^ a graphical representation of the H.F procedure. The thin lines 
stand for free polymer Green function Go(r, r), while the thick ones correspond to 
the self-avoiding polymer Green function G{f,T). The calculations performed in the 
present paper incorporate the contribution of the two upper diagrams. The correction 
associated with S*^^^(r2 — fi) will be very briefly discussed in the next section. 

2.4 The HF equation in momentum space. 

Because of the space translation invariance of eq. P5| ) considerable simplification is 
achieved by working within the plane wave basis where the operators G, Gq and S 
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are diagonal : 

{P2\G\pi) = {P2mT))^\pr) = 5\p,-n)G{p,,T) 

= P2) J d'^ r exp{ipi ■ r) G{r). (31) 

The Hartree-Fock integral equation can then be written under the rather compact 
form: 

(5(pi,r)-i = Go(pi,r)-i-S(pi,r) 

2 

= Go{pl,r)-' + J^|d'pGip,r)Vi\p^-p\). (32) 

In order to compute the extension vs force curve in the general case, one has to 
proceed as follows: 

i) First, we take the limit F = 0. The operators Ro{t) and (i?(r))^ are then 
rotation invariant and G{p,t) = g{p'^,T). The HF integral equation becomes one- 
dimensional with as kernel K{pi,p) = 2np'^ j\ dxV (y\Jpi — p"^ — 2 x pi pj . In the 

particular case of Yukawa or exponential potentials, K{pi,p) is easily computed an- 
alytically. Returning to Statistical Mechanics notations one sees on eq.(|27|) that 
gijP'jT) is nothing but the Laplace transform z{p'^,t) of Z{p'^,N), introduced in eq. 

(O)- 

ii) The second step involves the analytic continuation p"^ —F"^ of z{p'^, r). Then, 
one has to perform the inverse Laplace transform in order to get the extension vs force 
curve C{F, N) from eq . (|I8D . As we shall see, for the particular case of the 6 function 
potential, ({F, N) is dominated in the large N limit by the contribution of the pole 
of g{—F'^, t) = z{—F'^, r), occurring at the largest algebraic value of r. 

3 The HF Method for a Short Range Repulsive Monomer-Monomer 
Potential. 

From now on, we are going to use b as the unit of length and as before /c^T as the 
unit of energy. The short range limit of the Yukawa potential g"^ V{r) = |^ exp — ^, 
is obtained by taking a — 0, keeping g^ finite. The result is just the (5-function 
potential : V{r) = g"^ a? 5'^{f). The HF proper self-energy reduces to: E(r2 — "H.) = 
-g^ G(0, r)(53(f2 - n). This implies: E(p, r) = -^^ ^3 ^(g, r). 

3.1 The reduction of the HF integral equation to a numerical equation. 
The Green function in momentum space is then given by : 

G(p, r) = g{p\ r) = {h,{p') + r + g' a' G(0, r))'' . (33) 
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By computing G{f, r) from the above equation by a Fourier integral and then taking 
the hmit r — > 0, one gets a self-consistent equation giving ^(O, r): 

^^(0'^) = ^^^ /io(p2)+/+ j2a3G'(0,T) ■ ^^^^ 

Prom now on, we shall take /lo(p^) = h^'^\p^) — \p'^ -\- ^P^- ( Remember that 
our length unit is h). With our choice of hoip^) the integral in the r.h.s. is clearly 
convergent, while it would have been linearly divergent in the case of the free Gaussian 
chain. We have found convenient to use instead of G(0, r) the following self-consistent 
function: 

//(T)=T + /a=^G'(0,T). (35) 
To solve exactly the HF equation we have introduced the algebraic function : 

/■oo ^2 TT f-1 + 2 K Vr + Vl - 4 k2 r) 
F{i^:r)^ = ^ , ^ = , (36) 

where we have made explicit for clarity the cutoff parameter = / = 1/30. 
In order to study the cross-over occurring in the vicinity of / = , the following 
approximation of F{k, t) is useful: 

FoM-^^ (^l-«:V^+^ + 0(«:v^)'j. (37) 

We are now ready to write down the self-consistent equation giving /x(r) : 

h{t) = t + ^F(«,6/x(t)), 

A = ^/al (38) 

TT 

Remembering that g{—p, r) — z{—p, r) = (/io(— /^) + //(r))"''" is the Laplace trans- 
form the searched- for quantity Z[—f^,N), one has to perform the following inverse 
Laplace transform: 

Z{-f,N) = —- / dreMNr) (/io(-f ) + //(r)) . 

It appears that the relevant singularities of g{—f^,T), in the complex half plane 
9(t) < 0, are lying along the negative real half axis. As a consequence, the integration 
contour along the imaginary r axis can be folded around the negative real half axis. 
This leads to the following integral formula : 

Z{-f,N) = --{ + )rfrexp(iVr)(V-/')+Mr))" . (39) 
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Figure 4: The partition function Z{—f'^,N) is obtained from the inverse Laplace transform of 
(^o(~/^) + ^ ■ The result is dominated by the pole residues of the integral. Since /io(~/^) < Oi 

we have to look for the positive solutions of /^(r) = t + {K,6fi{T)). The continuous line gives 
the self consistent function /i(T) for A = 0.7. The dashed line is associated with the approximate 
function fJ-oir), obtained with the simplified function Fq{k,t), valid for small r up to corrections 
of order (55)^^^. The figure shows that this approximation is adequate in the region below the 
horizontal dotted line, which corresponds to the low force interval < / < 0.5 . 



3.2 Analytical computations of the extension vs force curves in the low force regime. 

Let us consider first the low force regime < / < 0.5 i.e < — /io(^/^) ^ 0.0413. 
We are clearly interested by the value of fi^r) ~ — ^o(~/^) < 0.0413. The use of 
the approximation 6;u(r)) ~ Fo(k, 6/i(r)) is justified here by the small value 
expansion parameter : K^Qjiij) < 0.091. The corresponding approximate function 
/io(T) differs from the exact one by less than 0.5% if < / < 0.5 , as it is seen clearly 

on FIG. ^. Let us look first for the value r = Tc = — — 2 — syxch that /i(rc) = 0. It is 
convenient to define the translated variable a = t — r^.. To get an explicit expression 
of [loir), what we have done first is to solve the self-consistent equation for r > 0. In 
this region , ^io{t) > , so that the function Fo(k, 6/io(T)) is regular. We perform next 
an analytic continuation toward the point t = Tc where f^o^r) vanishes. According to 
eq. ( refinvlap), we have to follow paths running just above ( and below ) the half 
negative half axis. This insures that y/JjF = /i in the vicinity of r = Tc. We arrive 
then to the following explicit formula, valid in the region of physical interest: 
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240a + 7rA ( 5 tt A - 6 a/30 a + ,/2400 (X + 5 tt A ( 5 vr A - 12 V30 (X 



24^2 144 f-40 + VSOvtAW^ 

Some remarks are in order : first, the above expression has for small value of r — Tc, 
the parabohc shape exhibited on FIG.^, second a branch point occurs for t — Tc = 
Tb = _4go+i2^7r A " nianifests its presence on the curves of FIG.§ by a spike. (Below 
the branch point it is actually 3fJ(/io(T)) which is plotted.) Such a singularity does not 
seem to show up in the range of interest for the exact function fi{T) and it is likely 
to be an artefact of the approximation F{K,Qfi{T)) ~ Fo(fi:, 6/i(r)). Nevertheless, 
we have computed explicitly its contribution to the inverse Laplace transform and it 
was found to be negligible compared to the dominant poles contributions. Finally 
it is of interest to note that the coefficient of the series expansions with respect to 
a are divergent in the weak coupling limit A ^ 0. Furthermore, we have checked 
the identity of series expansion in power of A obtained by two different ways: one, 
by a direct expansion of the above exact expression of fJ.oi'r), two, by solving the 
self-consistent equation by a perturbation method. The coefficients of the various 
powers of A have a branch cut at r = which disappear when the exact summation 
is performed. The above remarks provides an illustration, on a simple case, of the 
importance of non-perturbative effects in the self-avoiding polymers problem. 

The two poles ti 2(/, A) of (/io(— /^) + l^oiT)) ^ are easily computed: 

ri,2(/,A) = T,± +7]{f) (1 




vif) = f\n-^P- (41) 



We have in hands all the information needed to compute the partition function 
Z{—p,N) given by eq. (0). Ignoring the branch point contribution, we are left 
with the residue contributions relative to the two poles at r = ti 2(/, A): 

Z{-f\ N)=j: exp (AT n{f, A)) ( ^-^^MIM] " . (42) 



i=l 



The relative extension ({f,N) = {z{N) )/N is then given by the logarithmic deriva- 
tive: 

C(/,iV) = -^^log(^(-f,iV)). (43) 
18 
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Figure 5: On this figure we have displayed a set of extension vs force curves in the low energy 
regime, for typical values of the monomer number A'^. The bottom thick dashed curve represents 
the FJC extension vs force curve with no self-avoidance. The most remarkable feature is the steady 
increase with N of the slopes at the origin (/ = 0). We note also that in the limit of very large N, 
the extension jumps to a finite value at very small force. This suggests the presence of a cross-over 
phenomena near f — 0. This could be associated with the partition function Laplace transform 
poles crossing at the origin. 

Using the explicit formulas for //o(t) and rj(/, A) given in eq. ( ]iO| ) and (^T]), ({f, N) 
can be computed analytically. The corresponding expression is rather involved and 
not very illuminating. A better illustration of what is going on in the low force regime 
is obtained by looking at FIG. |] where relative extension vs force curves ({f,N) = 
{z{N) )/N are displayed for increasing values of the total chain monomer number N. 



3.3 Comparison of HF extension slopes at the origin (f=0) with the results derived 
from RG arguments. 

The most salient feature of these curves lies in the fact that the extension vs force 
curve slope, s(A^, A) = ^^j^|j=0; increases with while it stays constant in the 
case of the FJC model. This kind of behaviour is found under a less pronounced 
form in the Edwards model (Gaussian chain with short range self-avoidance), using 
renormalization group arguments (RG). We would like to stress that a quantitative 
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Figure 6: In this figure we present the results of a simple analysis of the variation of the slope 
of the extension vs force curve s{N,X) at f — with the monomer number N. We have plotted 
the characteristic quantity a{N, A) defined as N times the logarithmic derivative of the slope at the 
origin with respect to N. If a{N, A) stays constant, the slope behave as N" or, in other words, it 
obeys a scaling law. RG arguments applied to the Edwards model gives a = 2iy — 1, where v ~ 0.6 
is the Flory coefficient. Our model exhibits a moderate but significant variation of a{N, A) with N. 
Moreover, a{N,X) is not "universal ", in the sense that it is roughly proportional to the coupling 
constant A. We would like to stress that this difference was to be expected. As we have pointed out 
previously, our model can be viewed as a fixed cutoff ( = 30/5^) version of the Edwards model 
and, as a consequence, RG arguments are not valid. 



comparison of the two approaches is not a simple problem. We have noted previously 
that the model used in the present paper can be viewed as a version of the Edwards 
model, endowed with a a fixed cutoff = -p. It is then clear that under such 
circumstances RG arguments are not appropriate. 

For exhibiting the difference of behavior in the vicinity of / = 0, it is convenient 
to introduce the following quantity : a{N,\) = N log {s{N, \)). Its physical in- 
terpretation becomes more apparent by noting that if a{N,X) is constant, then the 
slope s{N,X) is proportional to A^". This corresponds to the scaling prediction |]16 



for the Edwards model with a = 2z/ — 1, where z/ ~ | is the Flory coefficient, which 



governs the variation of the root mean square radius of the Edwards chain in absence 
of stretching force : { r^(A^) ) oc N'^. Curves giving the variation of a{N,\) are 
displayed on FIG.^. In the case of the present HF computations, there is, for fixed A, 
a moderate but significant increase of a{N, A) with A^. Moreover, for fixed A^, when 
A goes from 0.2 to 0.7, a{N, A) is multiplied by a factor varying from 3 to 2 when N 
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increases from 50 to 500. 

These features are likely to be associated with the fact that we are dealing with 
a fixed cutoff field theory problem but they could also be partly due to an artefact of 
the Hartree-Fock approximation. In order to check this latter point, one should try 
to estimate the corrections to the HF approximation. In the particular case of the 
short range repulsive potential g'^V{r) = a^g^S{r), the proper self-energy correction 

(r2 — 'f^i) is given by the rather simple expression: 

- ri) = A' (G(f2 - n)f (44) 

The self-consistent integral equation cannot be reduced anymore to the solving of a 
numerical equation. The resulting non-linear integral equation can be solved by an 
iteration procedure starting from the HF approximation. Such a calculation, which 
looks feasible, is clearly of considerable interest but it falls outside the scope of the 
present exploratory paper. 

4 Comparison of the Hartree-Fock Force vs Extension Curves Compu- 
tations with the Results of Monte-Carlo Simulations 

In this last section we would like to compare the results of the Hartree-Fock method 
applied to the self-avoiding effects for a continuous version of the FJC model with 
those obtained by a Monte-Carlo simulation of a self-avoiding freely joining chain. 

We begin by saying some words about the way we have obtained the theoretical 
force vs extension curve displayed on FIG.^. In the range: < / < 0.25 we have used 
the analytic formulas of the previous section, which are valid in the low force regime. 
In the medium force interval : 0.2 < / < 1.5 the computation was performed with the 
exact solution /i('r) of the self-consistent equation (^). A considerable simplification 
was achieved by noting that, if / > 0.2, the pole residue at r = ri(/. A) > r2(/. A) 
gives the dominant contribution, by more than two orders of magnitude. 

Above / = 1.5 we are leaving the domain of strict validity of our computation. We 
have performed an extrapolation to cover the interval 1.5 < / < 3.0. It corresponds 
to the dashed part of the theoretical curve of FIG.|^ We do not extrapolate directly 
the final result but rather the quantity g{—piT)~^ = ho{—f'^) + /^(t). For the first 
term, we set ho{—f^) = /iq(— /^), instead of keeping h^^—p) which would lead to 
an inadequate force vs extension curve for / > 2, as it appears clearly on FIG.|I|. 
Concerning yu('r), we note two things: first, the replacement h^lp"^) ^o(p^) will 
only affect the high p tail of the convergent integral giving the function F{k,, t) (see 
FIG.0), second, the non corrected /i(T) varies linearly in the interval of interest: 
0.35 < h{t) < 1.2. It look then reasonable to perform our extrapolation by keeping 
the self-consistent function /i(r) given by eq. (|38|) . 
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Figure 7: This figure presents a comparison of the force vs extension curves for N = 100 obtained, 
on one hand, from the HF analytic computations of the present paper, on the other hand, from MC 
simulations performed by the ENS experimental group Our analytic treatment, represented by the 
continuous curve, is strictly valid in the low to medium force range < / < 1.5. Nevertheless, we 
have performed a reasonably safe extrapolation - the dotted branch of the curve- in order to cover 
the higher force range 1.5 </ < 3. The choice of the value of our single parameter: A = 0.7 does 
not come out from a least square fit but from a trial and error procedure. In the MC simulations 
self-avoidance was implemented by keeping only the configurations such that the distances between 
any two non -adjacent monomers are larger than a fixed length d. In the actual simulations b was 
taken equal to the fixed distance b between adjacent polymers. The fairly good agreement between 
these two different approaches is quite remarkable and it may suggest the shape of the force curves 
is not very sensitive to the range of the repulsive potential if d < b. 



The M.C simulation results displayed as big dots on FIG.0 have been obtained 
by the authors of reference |T^, with a chain of = 100 monomers. Two adjacent 



monomers are freely connected by a segment of length b. The repulsive potential 
between two non-adjacent monomers is taken to be of the square wall type: V{\ri — 
^2!) = Vo if |ri — r2|) d and otherwise. In the simulation, it is tacitly assumed 
that Vq ^ ksT . In practice, this means that the only accepted configurations 
are such that the distances between any pair of two non- adjacent monomers are 
> d. The short range potential used in our calculation corresponds to a different 
hmit : Vo/kBT ^ 00 , d ^ 0, with ^d^Vo/kBT going to a finite limit -fb^. The 
limiting potential V{\fi — r2|) reduces to the zero range potential jb^Slfi — r2). The 
coupling constant A used in our computation is to be identified with ^76'^. Despite 
the difference of the two kinds of limit, we see on FIG.0 that the M.C. simulation 
results with d = b are fairly well fitted by our HF curve, if we choose A = 0.7. 
( We work in a system where b and kgT are taken respectively as units of length 
and energy). In physical terms, this may suggest that the shape of the force vs 
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extension curve is rather insensitive to the potential range d, as long as it is equal 
or smaller than the monomer length b. This important point could be confirmed in 
two ways: one by MC simulations with d/b < 1, second by solving numerically the 
HF integral equation (^), involving finite range Yukawa or exponential potentials. 
If this last numerical program can be achieved successfully, it will also open the road 
to HF method applications to biological polymer like ssDNA, using more realistic 
monomer-monomer potentials. 



Conclusion 

The exploratory investigation of the HF method, as a tool for the study of fiexible self- 
avoiding polymers, appears encouraging enough to justify an extension of the present 
work in several directions. The most urgent task is perhaps to get an estimate of 
the "non planar" diagrams contributions to the self-consistent equation, along the 
lines suggested at the end of section 3.3. Despite the good agreement with the MC 
simulations, it is far from obvious that the main physical results of the present work 
stay robust vis-a-vis such corrections. 

The second road to be explored is the numerical solving of the self-consistent 
equation ( P^j) with more realistic monomer-monomer potential V{r), say a Yukawa 
potential, rather than the short range limit used in the present paper. The numer- 
ical work load will not be sensibly increased if one takes for V{r), a superposition 
of Yukawa and exponential potentials, adjusted in such a way that the monomer- 
monomer interactions are attractive for distance about the hydrogen bound length 
and become repulsive when one reaches distances in the range of the Debye length. 
The MC simulations and the HF model force vs extension curves deviate from the 



data relative to inhibited hydrogen bonding |T0[, in the low force regime. In par- 
ticular the experimental curves suggest that the molecular chain extension vanishes 
for a small but finite value of the stretching force, about 0.05pN. This could be an 
indication that hairpin structures are still present at very low force, say below O.lpN. 
It will be of interest to see, if by adjusting the strength of the attractive monomer- 
monomer interactions, our HF model is able to reproduce the data in the whole range 
of force. 
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